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O ABSTRACT 
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We show that main-sequence stars in dense stellar cusps around massive black holes 
, are likely to rotate at a significant fraction of the centrifugal breakup velocity due to 

j. \ spin-up by hyperbolic tidal encounters. We use realistic stellar structure models to 

t-h | calculate analytically the tidal spin-up in soft encounters, and extrapolate these results 

to close and penetrating collisions using smoothed particle hydrodynamics simulations. 

■ We find that the spin-up falls off only slowly with distance from the black hole because 

■ the increased tidal coupling in slower collisions at larger distances compensates for 
T ■ 

^sj . the decrease in the stellar density. We apply our results to the stars near the massive 

black hole in the Galactic Center. Over their lifetime, ~ 1 Mq main sequence stars in 
the inner 0.3 pc of the Galactic Center are spun-up on average to ~ 10%-30% of the 
centrifugal breakup limit. Such rotation is ~ 20-60 times higher than is usual for such 

■ stars and may affect their subsequent evolution and their observed properties. 

6 ■ 

Subject headings: galaxies: nuclei — Galaxy: center — Galaxy: kinematics and 
dynamics — stars: rotation — stars: kinematics 
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■ 1. Introduction 

It is now widely accepted that super-massive black holes (BHs) exist in many, if not all 
galactic centers (Magorrian et al. 1998j ). Dynamical models of the evolution of such systems 



generically predict the formation of a dense stellar cusp near the BH (e.g. Bahcall & Wolf 1977 



Young 1980| ). Stars moving rapidly in the dense stellar cusp near the BH will suffer numerous 



high-velocity close tidal encounters over their lifetimes. Although such encounters transfer some 
energy and angular momentum from the hyperbolic orbit to the colliding stars, they rarely remove 
enough energy for tidal capture. This is in marked contrast to the situation in the high density 
cores of globular clusters, where the colliding stars are on nearly zero-energy orbits and close 
collisions lead to the formation of tight binaries. The effects of hyperbolic encounters on the stars 
are mostly transient. The dynamical and thermal relaxation timescales are very short compared to 
the stellar lifetime, and thus apart from some mass-loss in very close collisions, the star is largely 
unaffected. It is however more difficult for the star to shed the excess angular momentum since 
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magnetic breaking operates on timescales of the order of the stellar lifetime (Gray |1992 ). High 
rotation is therefore the longest lasting dynamical after-effect of a close encounter. 

The possibility that stars in BH cusps are rapid rotators may have interesting implications 
for their evolution and the interpretation of their observed properties. The effects of rotation 
and rotationally induced mixing in a main-sequence (MS) star on its subsequent evolution has 
been studied by Sweigart (1997) and VandenBerg, Larson & De Propris ( |1998 ) in the context 



of globular clusters, without specifying the origin of the rotation, and by Sills, Pinsonneault & 
Terndrup ( [2000| ) in the context of young open clusters. Rotationally induced mixing may reveal 
itself in the spectral line ratios, and rotation may be directly observed in the spectral line profiles. 
Detection of such signatures in the spectra of the observed giants can provide additional evidence 
for the existence of an underlying cusp of MS stars, which at present cannot be directly observed. 

The goal of this study is to estimate the magnitude of the tidal spin-up, with particular 
emphasis on the Galactic Center (GC). Present-day observations can already resolve individual 
giant stars very close to the BH in the GC (Genzel et al. |1997 ; Ghez 1998) and high resolution 



infrared spectroscopy is possible for the brighter giants (Carr, Sellgren & Balachandran 199E; 



Ramirez et al. 2000). The depletion of luminous giants in the inner 2" (~0.1 pc) around the 



supermassive BH in the Galactic Center (GC) was interpreted by Alexander ( |1999j ) as evidence for 
collisional destruction in an extreme density of a sharp stellar cusp. It is inevitable that where the 
stellar density is high enough to destroy giants, smaller stars that escape destruction will suffer 
very close collisions. The inner GC is therefore a promising environment for studying the spin-up 
effect. 

The paper is organized as follows. In Sec. ^ we present the formalism for spin-up in the 
linear regime of soft hyperbolic encounters. In Sec. ^ we discuss results from SPH simulations 
of non-linear close encounters and incorporate them in our spin-up calculations. In Sec. ||] we 
calculate the spin-up of stars in the inner parsec of the GC. We discuss and summarize our results 
in Sec. ||. An appendix describes the analytic calculation of the tidal coupling constants for 
hyperbolic encounters using realistic stellar structure models of a MS dwarf and of a giant. 



2. Spin-up by soft hyperbolic tidal encounters 

2.1. Stochastic spin-up in the linear regime 

We begin by considering soft (distant) encounters where the tidal deformations are small 
enough to be treated as linear perturbations. We consider the effect of the tides raised by an 
impactor star of mass m on a target star of mass M and radius R as the impactor follows an 
unbound orbit with a peri-separation r p from the target star. We will use the tilde symbol to 
denote mass in terms of M, distances in terms of R, time in terms of y/R 3 /GM , velocity in terms 
of the Keplerian velocity y^GM/R, energy in terms of GM 2 /R, angular momentum in terms of 
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V GM 3 R and moment of inertia in terms of MR 2 . In these units ft = 1 is the centrifugal break-up 
angular frequency. A star with $7 > 1 will shed mass from its equator due to the centrifugal force. 

We describe the results in the reduced mass system where the target star is at the origin. 
The angular momentum AL that is transferred from the orbit to the target star is related to the 
deposited tidal energy AE by (Goldreich & Nicholson 198S ; Kumar & Quataert 1998| ) 



AE = Am p , (1) 

where f2 p is the angular velocity at periastron. The energy invested in raising the tides is given in 
the linear regime by (Press & Teukolsky |1977| ) 



r P 1=2 r P 



where the tidal coupling coefficients 2} depend on the structure of the star and the eccentricity of 
the orbit and are functions of the parameter ij, 



r) = JjL r . (3) 

For rigid body rotationF], the change in the star's angular momentum is related to the change in 
the angular velocity Afi by the star's moment of inertia /, 

AL = I A^l . (4) 

In the case of differential rotation, Af2 is defined by Eq. || and is the effective angular velocity. 
The spin-up of the target star in a single tidal encounter can therefore be expressed as 

A S = J^£?#, (5) 

llL p'p 1=2 P 

where it is assumed that the star maintains its original moment of inertia. The periastron angular 
velocity is related to the relative velocity at infinity, Vqo by 

~ 2 t& 2(1 +m) 
P P 

where the second term expresses the enhancement due to gravitational focusing. 

Many encounters at random orientations lead to a "random walk" buildup of the stellar spin. 
Over the stellar lifetime T, the rms change in the stellar angular velocity is given by adding in 



lr The timescale for angular momentum re-distribution due to convective transport in a late type giant is ~ 1 



(Zahn 1989) and so rigid rotation is achieved on a timescale similar to that of the collision itself. The timescale for 
angular momentum re-distribution in radiative MS stars is not well known, although it is likely to be shorter than 
the stellar lifetime. 
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quadrature the contributions from collisions with different values of the peri-separation r p and of 
the orbital energy E a , weighted by the differential collision rate d 2 q/dr p dE a , 

The differential rate is calculated with the approximation that the relative velocity can be 
described by the Maxwellian distribution function (DF) with a mass-independent ID velocity 



dispersion o\ = 2cr (Sec. 4.1) 



dr p dE Q H 2 4 \ p I 

where n is the space density of the impactors. This estimate involves the approximations that 
the angular momentum transfer does not depend on whether the star is rotating in a prograde 
or retrograde sense with respect to the orbit, and it assumes that the mass-loss and the change 
in the star's structure and moment of inertia can be neglected. Because Ail 2 ~ r~ 9 whereas 
d 2 q/dr p dE a ~ r p , the rms 5£1 is dominated by the collisions with the smallest r p . When the stellar 
population includes a spectrum of masses, / f~dfh = 1, the average spin-up is obtained by adding 
the weighted contributions in quadrature 

8ti = ( [ dmf~Atl 2 ) 1/2 . (9) 



2.2. Linear tidal coupling coefficients 

The orbital energy E Q is related to Voo and the eccentricity of the orbit e by 

l™o l_e-l 

where Jl = m/(l + m) is the reduced mass. The eccentricity that corresponds to the mean orbital 
energy (Eq. ^2] below) is 

(e) = 6r p a 2 + 1 . (11) 

Since e can reach very high values when a 3> 1, as is the case for giants (low Keplerian velocity) 
close to the BH, it is necessary to extend the standard (e = 1) tidal interaction formalism of Press 
& Teukolsky ( 1977| ) to hyperbolic orbits. This is described in Appendix |A|. 



In this work we consider two types of target stars: A MS dwarf and a red giant, whose 
properties are summarized in table |]. The detailed stellar structures that we use to calculate the 
tidal coupling coefficients for these stars are based on the solar model of Christensen-Dalsgaard et 
al. ( |1996| ) and a model for the red giant aUMa (Guenther et al. 20001) . Figure |l] shows the run 



of T\ with r] for the two stars and for an ideal gas n = 1.5 polytrope, for parabolic (e = 1) and 
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hyperbolic (e = 10) orbits. We find that the tidal coupling in an ideal gas polytrope is stronger 
than in either of the more realistic models. We also find a general trend for the T\ coefficients to 
reach their maxima at larger values of r\ (larger r p for a given in) with increasing e (and increasing 
periastron velocity v p ). This reflects the fact that the coupling is strongest when £l p = v p /r p 
equals the lowest frequency stellar pulsation mode. 



3. Spin-up by strong hyperbolic collisions 

3.1. SPH simulations 

In order to extend the linear treatment of the soft encounters to strong (close and penetrating) 
encounters, we simulated such encounters with the Smoothed Particle Hydrodynamics (SPH) 
technique (Lucy |1977 ; Gingold & Monaghan 1977). In view of the difficulty of simulating all the 



aspects of a real stellar collision, and in view of the many uncertainties in the details of the stellar 
structure, the purpose of these simulations is not to calculate AO, precisely for specific collisions, 
but rather to gain qualitative insight about angular momentum transfer in strong encounters, 
which can then be incorporated in our semi-analytic calculations by simple approximations. The 
SPH code we use calculates the gravitational force by straight N 2 operations, and is therefore 
limited to relatively low resolution simulations (typically N = 2048 particles). The code integrates 
in time the entropy equation (Hernquist |1993|) ; conserves particle momenta identically; uses the 



artificial viscosity prescription given by Hernquist & Katz ( 1989 , Eqs. 2.22, 2.23, 2.37) and the 



time-step criteria of Katz, Weinberg & Hernquist ( 1996| ). The amount of stellar mass- loss in the 



collisions is estimated by the enthalpy criterion of Rasio & Shapiro ( 1991 ) 



We verified the code by constructing stable n = 1.5 polytrope configurations; by confirming 
that the results converge as the number of particles is increased; and by reproducing qualitatively 
the spin-up and mass-loss obtained in the SPH simulations of Davies, Benz & Hills ( |1991| ), who 
used a much more realistic stellar structure model. Of direct relevance is the fact that our SPH 
code reproduces the results of the linear theory at the soft collision limit (r p ^ 2.5) and at smaller 
r p follows closely the SPH results obtained by Rasio & Shapiro ( |1991| ) (Fig. |2|) in SPH simulations 
with 10 4 particles. 



3.2. Beyond the linear regime 

3.2.1. Deep inelastic collisions 

At small peri-separations (small rj for fixed m), the sum in Eq. |5] converges slowly and the 
truncation of the 7} series at some order I = k could under-estimate AO. We make use of the 
fact that the ratio Ti + i(rj) /Ti{rj) is roughly constant (Fig. |l]) over the small-77 range of interest to 
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extrapolate the sum to high I by a geometric series 

fh 2 f k ~ 2 T T 

A n « CnL~~— £ ~2l + Jk^2 



L P' p \l=2 V 



(12) 



where the constant Cnl is a non-linear correction factor, which we calibrate by the SPH 
simulations as discussed below. Figure § compares the prediction of the linear theory with the 
numeric SPH results without correcting for non- linear effects (CVl = 1)- While the two agree at 
the limit of soft (distant) encounters, AE grows faster than the linear theory at close encounters. 
In this particular example the high-order correction is very small, and does not exceed 2% down 
to f p = 1.6. 

Figures. |3| and [| show snapshots from two SPH simulations of extremely non-linear collisions. 
Figure || shows the run of the dynamical properties of the target star with r p in fh = 1, E = 2 
collisions. Our analytic expressions for AL and AQ consistently underestimate the SPH results 
for a point mass impactor by a factor of ~ 3 down to r p ~ 1.5. Similarly, it under-estimates AE in 
parabolic collisions by a factor ~ 1.5 at r p = 1.8 (Fig. ^). We find that at smaller peri-separations 
where significant amounts of mass are ejected in the collisions, a large fraction of the angular 
momentum that is taken out of the orbit is deposited in the ejecta rather than in the target star. 
This truncates the AO, ~ r~ 5 divergence and limits the spin-up efficiency, so that the difference 
between the analytic calculations and simulation results decreases with r p until the two cross over. 
We approximate the decrease in spin-up efficiency by introducing a truncation peri-separation ro 
such that AQ(r p ) = Afi(ro) for r p < ro. 

We now discuss our choice of ro. Figure |] shows that AO, is much more suppressed than 
AL at small r p because the stellar moment of inertia increases as a consequence of the collision. 
The increase in / far exceeds the small change that is expected due to the stellar oblateness 
that develops in response to the rotation. The effect seen in the simulation is due to the heating 
and subsequent expansion of the star by the collision. Our SPH models do not include the 
radiative processes that are necessary for describing the later stages of cooling and contraction. 
It is therefore likely that the final value of AQ is not as strongly suppressed as is implied by 
Figure ||, but rather follows more closely the behavior of AL. A collision with an impactor of a size 
comparable to that of the target results in much more mass-loss than a collision with an impactor 
that is effectively a point mass (e.g. a stellar remnant on a MS star or a MS star on a giant). 
However, it is likely that the SPH results over-estimate the mass-loss since a n = 1.5 polytrope has 
a significantly lower binding energy than the realistic stellar structure models on which we base 
our analytic calculations (Table |l|). Figure || shows that the analytic estimate of AO, equals the 
value derived from the simulation at r p ~ 1.3 and r p ~ 0.9 for polytrope and point mass impactors, 
respectively. For AL the analytic estimate and the simulation results are equal at r p ~ 0.9 and 
r p ~ 0.8, respectively. Based on the arguments presented above, we adopt in this work the simple 
prescription that ro = 1.0 or the size of the impactor, whichever is larger, for all types of collisions. 
We note that the analytic calculations still under- estimate the SPH results by a factor of 2 at 



-7- 



ro, and by a factor of 3 at r p = 2. A difference of a similar magnitude is seen also in parabolic 
collisions (Fig. Q). We compensate for this discrepancy by setting the non- linear correction factor 
to Cnl = 2. A more extensive investigation of parameter space by SPH simulations with more 
realistic stellar structure models will be required to refine this prescription. 



3. 2. 2. Prompt disruption of tidally formed binaries 

Close hyperbolic encounters of stars in the low energy tail of the orbital energy distribution 
can lead to the formation of a bound system {E Q + AE Q < 0) with a large semi-major axis a. 
Although such encounters are very efficient in spinning- up stars, the subsequent evolution of the 
stellar rotation in binaries is very different from that due to stochastic encounters. We do not 
consider such cases in this work and they are not included in the average <50 (Eq. ^). Because 
of the very large stellar density and the proximity of the massive BH, it is necessary to check 
whether a newly formed binary can survive its first orbital period without being disrupted by a 
tidal interaction with either a third nearby star or the central BH. 

We take this into account by considering the orbit as bound only if the change in the orbital 
angular momentum, AL or b, due to the torque exerted on the system by the third mass (star 
or central BH) is smaller than the orbital angular momentum -L or b, 

a\ 3 ( fhz \ / o\ -1 / 2 



where d is the distance of 777,3 from the binary and where we used the relation -L 2 rb = fhjia(l — e 2 ). 
For disruption by a star, d ~ n -1 / 3 , fh^ ~ fh and the no-disruption criterion is 

4vrna 3 /I (l - e 2 ) < 1 . (14) 

The eccentricity and the semi-major axis are estimated from the unperturbed peri-separation and 
the bound orbit's energy 

H 1 — e 

In practice, we find that the contribution from disrupted binaries to <5f2 in the GC is negligible. 



3.2.3. Tidal and collisional destruction 

The precise criterion for tidal destruction of stars by hyperbolic encounters is not well 
known. In this work we adopt the simple criterion that tidal break-up occurs when Av, the 
change in the velocity of a mass element on the stellar surface, exceeds the escape velocity, 
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Av ~ 2m /r p x r p /v p > \[2. The tidal radius Ft is a function of the impactor mass and orbital 
energy and is given by the solution to the equation 

/ \ 1/3 

r t =[-^) m^. (16) 
\m + r t E Q J 

Our SPH simulations indicate that this is a conservative criterion. Figure [5] shows a sequence 
of snapshots from an encounter with r p = r t = 2.0, E Q = 1.0 between a m = 10 black hole 
(modeled as a point mass) and a MS dwarf (modeled as an ideal gas n = 1.5 polytrope). The 
star survives the collision after suffering ~ 25% mass-loss. We performed several spot-checks with 
different values of r p and E Q and verified that this prescription for Ft indeed roughly demarcates 
the boundary where the fractional massdoss increases to order unity, and that A$7 significantly 
exceeds the value predicted by our extrapolated linear formalism for such deep collisions. 

When Ff < 1 (point mass impactor) or Ft < 2 (impactor of same size as target) it is necessary 
to consider the possibility of collisional destruction. Our SPH simulations indicate that the n = 1.5 
polytrope stellar models can survive very deep collisions (r p ~ 0.5, cf Figs ||, |j), albeit with a 
significant massdoss. One concern when considering penetrating collisions by compact remnants 
is the energy release by nuclear reactions near the surface of the impactor. Rbzyczka et al. ( |198S| ) 
and Ruffert & Midler ( 1990| ) find that nuclear reactions probably do not play a significant role 



even in parabolic head-on collisions between low-mass stars and white dwarfs. We will assume 
that stars can survive collisions down to r p = 0.5 even when the impactor is a compact object. 



3.2-4- The survival probability 

In an environment that is dense enough for efficient tidal spin-up there is also a non-negligible 
probability for destructive head-on collisions. Our estimate of Afi implicitly assumed that the 
target star survives its full life span T. We now estimate the survival probability of stars against a 
close collision with some peri-separation F c in order to check what fraction of the star can survive 
long enough to acquire significant rotation. For simplicity, we omit gravitational focusing, which 
is negligible for long-lived MS stars very close to a massive BH (a ^> 1). 

The collisions are a random process in both r p and in the time of the periastron passage. The 
rate for collisions with peri-separation in the range r p to r p + dr p is 

dq = nVoo2T:r p dr p , (17) 

where in a Maxwellian DF with ID velocity dispersion a the rms Foo is \^6a. The probability 
density function p c (pdf) for a star to have over a time interval T c a closest encounter with 
peri-separation r p in the range [F c , F c + dr c ] is given by 

p c dr c = exp ( --nnvoofcrl ) 2nnv OQ fr c dr c , (18) 
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where the first term is the Poisson probability for avoiding an encounter with r p < r c and the 
second is the probability for having at least one in the required range. The fraction of target stars, 

f c , that avoid a collision with peri-separations r p < r c (the survival probability) can be written in 

-1/2 



terms of (r c ) = [AnVoaTA as 



POO 

fc = L Pcdr p = exp 

J r r 



IT 
1 



(rc) 



2 ) • (19) 



4. Tidal spin-up in the inner Galactic Center 

Up to this point our treatment of the tidal spin-up effect was general. We now turn our 
attention to the specific case of the GC. 



4.1. The stellar velocity distribution in the GC 

The effects of the stellar collisions depend critically on the relative velocity of the two colliding 
stars and their mass ratio. It is therefore important to understand the mass dependence of the 
velocity distribution. We are interested in particular in the case where the stellar system around 
the black hole has undergone two-body relaxation, as appears to be the situation in the GC 
(Alexander [1999| ). 

Let / m (e) be the DF of stars of mass m (mi < m < 777-2) as function of specific energy in a 
spherical stellar system whose potential is dominated by a central BH, where e = ^ — v 2 /2 > 0, 



\& = GM,/r, and M. is the mass of the BH. Bahcall and Wolf (|1977|) have shown that 
when such a system undergoes two-body relaxation, the DF has the following properties: 
p m = dlnf m /dhie ~ const, to a good approximation; p m2 ~ 1/4; and Pm/Pm 2 — Tnjm-i so that 
<s p m <j 1/4. In this case f m (e) oc e Pm and the velocity dispersion a m = (v 2 ) /3 is 

n / 1 \ GM. , . 

\Pm + 5/2J r 

The velocity dispersion is almost independent of the stellar mass. The relative change in the 
value of a m over the full mass range is only ~ 10% independently of the ratio 7772/7771, in marked 
contrast with the wide spread of velocities expected in the case of equipartition where a m oc m -1 . 
The reason why the relaxed system does not reach equipartition can be understood by considering 
the fate of a massive star that is momentarily on a circular orbit. The orbital radius depends only 
on the specific energy. Equipartition works to equate the kinetic energy per star, thereby always 
reducing the specific energy of the massive stars and causing them to sink to ever lower orbits. 
This is analogous to the equipartition instability discussed by Spitzer ( |196S| ) in the context of a 
stellar cluster without a central BH, but where the fraction of mass in the massive stars is large 
enough to create a centrally concentrated sub-system. 
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The resulting distribution of the relative velocity i>2 = v a — Vb between two stars of masses 
m i < m a < nib < rri2 is given by 

/w _ j_ £(i±^£(i±^j / A f 1 _ ^r- a _ (-z-jrp , (21) 

2 ' 8^*3/2 r(i+j> m jr(i+ Pm j J ^ 2^ ^ 2 / 

where u = t) a /^ 1//2 and = v-2,1^^ 2 and the integration is over the region u < \/2 and 
\w — u\ < \pl. This distribution is not very different from a Maxwellian DF with the same ID 
velocity dispersion a. This can be seen in Fig. 0, which compares the Maxwellian DF with the 



case where both p ma = p nib = (low mass stars) and Eq. 21 simplifies to 



f(v 2 )dv 2 = w 2 (32V2 - 24w + w 3 ) dw . 



24vr 

In particular, the two DFs have the same mean orbital energy 

(E ) = 3fia 2 . (22) 

We conclude that the velocity field in a relaxed stellar system very near a BH is well 
approximated by a Maxwellian velocity distribution where the ID velocity dispersion is 
independent of the stellar mass. 

4.2. The stellar population in the GC 

Stellar population synthesis models of the observed luminosity function averaged over the 



inner few parsecs of the GC (Alexander & Sternberg 1999|) indicate that it is well described by a 



continuous star forming population with a Miller-Scalo initial mass function (IMF) (Miller & Scalo 



1979) with masses in the range 0.1-125 M@. The mean impactor mass (live stars and remnants) in 
this model is (m) = 0.5 M@. In an isolated system the present day mass function (PMF) preserves 
the IMF distribution for low mass stars that are longer-lived than the system, but falls more 
rapidly for the shorter-lived massive stars, so that over time the low-mass stars accumulate and 
take an ever larger fraction of the total stellar mass. However, the timescale for mass segregation 
in the GC, which is of the same order as the relaxation timescale, is only 3Gyr (e.g. Alexander 



1999). It is therefore reasonable to assume that the mass fraction of the very low-mass stars in the 



inner GC is significantly lower than implied by continuous star formation. 

In view of the uncertainties in the low-mass end of the PMF, we do not attempt to construct 
a detailed mass function for the impactors in the innermost GC. Instead, we base it on the 
stellar synthesis model (table ||) and take account of the mass-segregation by assuming a Salpeter 
power-law PMF with a low-mass cutoff close to the mean impactor mass. The model parameters 
are listed in table ||[ The mean impactor mass in this model is (m) = 0.9 Mq. We carried out 
spot-checks to verify that the exact values of the mass ranges and power-law indices do not affect 
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the final results significantly. The one important assumption in this model is that the mass in the 
inner GC is not dominated by very low-mass stars, which are inefficient in raising tides. More 
detailed modeling of the dynamical evolution of the inner GC will be required to verify this. The 
results for the tidal spin-up that are presented below can easily be scaled to other mass fraction 
ratios through Eqs. || and |9[ 



We represent the mass distribution in the innermost GC by a stellar cusp of the form (Sec. 4.1 
here; Alexander 1999| ) 



-(3/2+p m ) 

10 77^ M pc- 3 , (23) 



,0.4pc 

and the ID velocity dispersion by Eq. [2(]. We set p m = to represent the typical low-mass 
impactors and assume M. = 2.6 x 10 6 M (Genzel et al. |l997|) . 



4.3. Results 

We calculated the spin-up in the GC following the procedure outlined in the previous sections. 
To summarize, the calculation proceeded by a triple numeric integration. First, integration 
over the impactor mass function (table ||). Second, integration over r p from a suitably large 
distance down to the larger of Ft (Eq. |l~6| ) and r p = 0.5 where for r p < 7q = 1, AQ was held 
fixed to account for mass-loss. Third, integration over E Q from a suitable large value down to 0. 
Collisions that resulted in bound orbits (Eq. |l^) were omitted from the sum. For each point in the 
integration, Af2 was calculated (Eq. |l^) with a non-linear correction factor Cnl = 2 and summed 
in quadrature, weighted by the differential rate (Eq. We assumed that long after the collision 
/ recovers its initial value and that fh is unchanged. 

Figure |7| shows the run of SO, with distance from the BH in the GC for the model MS 
target star with T = 10 Gyr, as well as the separate contributions from collisions with white 
dwarfs (WDs), neutron stars (NSs) and stellar mass BHs. Figure |7] also shows that the survival 
probability is almost unity as close as 0.02 pc to the BH, and so the fact that it was not taken into 
account explicitly in the estimate of 5£l does not introduce a serious error. The stellar rotation 
falls only slowly with distance from the black hole and is at the level of 5£l ~ 0.1-0.3 in the inner 
0.3 pc. Most of the effect comes from collisions with MS star and WDs. Field MS stars later 
than ~F5 are very slow rotators (Gray |1992|) . For ~1M Q stars these values of 5Q correspond to 
rotational velocities 20-60 times higher than normal (Table ||). 

We carried out similar calculations for giant stars. We find that the spin- up is much smaller, 
50 ~ 0.01-0.02, because of the much higher a and the shorter lifetime in the giant phase. This 
is similar to the rotational velocity of field giants later than ~ G3, and thus tidal spin- up may 
increase their rotation to double the normal value. 
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5. Discussion and summary 

Dense stellar cusps around massive BHs are environments where stellar collisions are frequent 
and energetic. Hyperbolic head-on collisions may destroy the colliding stars, but are rare. For 
every head-on collision there are many more close tidal encounters and grazing collisions. The 
cumulative effects of such encounters are more subtle. Energy and angular momentum are 
transfered from the orbit to the stars, during the fly-by the stellar structure may be significantly 
disturbed by the tidal forces and some mass may be lost. However, the stars survive the collision, 
and because of the high initial orbital energy, they rarely form a bound binary system. This is in 
contrast to the situation in dense globular clusters. The stellar dynamical and thermal relaxation 
timescales are very short relative to the stellar lifetime and the star can radiate the excess energy 
quickly. It is harder to shed the excess angular momentum, since magnetic breaking operates on 
timescales of order of the stellar lifetime. Thus the direct long term effects of the collision are 
increased rotation, and possibly some mass-loss and some mixing of the stellar envelope. 

In this study we calculated the magnitude of the rotation that is built up in a random walk 
fashion as the star undergoes multiple hyperbolic tidal encounters. This is of interest because 
of the effects high rotation may have on the star's evolution and on its observed properties and 
because, unlike destructive collisions, the spin-up affects the entire stellar population and extends 
over a much larger volume of the galactic nucleus. Many, if not most galaxies have a super-massive 
BH in their nucleus, and so high stellar rotation in galactic nuclei may be common. 

Our approach to the problem was to use detailed stellar structure models to calculate, for 
the first time, the tidal coupling constants for arbitrary hyperbolic orbits in the linear regime 
of soft encounters. Because the tidal energy and angular momentum fall off as a high power of 
the peri-separation, it was necessary to extend the calculation to the strongly non-linear regime 
of grazing and penetrating collisions. We carried out a suite of SPH simulations to study the 
qualitative behavior of such collisions and then conservatively extrapolated the exact linear 
calculations to the non- linear regime by several simple prescriptions. The simulations indicated 
the following. (1) At peri-separations closer than about twice the target star's radius, the linear 
results are smaller than the actual spin-up by a factor of at least two. (2) Once mass-loss becomes 
significant, the ejecta carries away a large fraction of the angular momentum that is extracted 
from the orbit, and therefore the spin-up saturates at its level just before the onset of mass-loss. 
(3) Stars can survive deep collisions down to a peri-separation of half the target star's radius. (4) 
A simple tidal disruption criterion for hyperbolic encounters can roughly indicate the point where 
the fractional mass-loss increases to order unity. 

We calculated the spin-up of stars in the inner parsec of the Galactic Center, which is of 
special interest because of the high quality and wealth of details of the observed stellar data. We 
find that over 10 Gyr, ~ 1 Mq MS stars in the inner 0.3 pc are stochastically spun-up to 10-30% 
of the centrifugal break-up velocity (20-60 times higher than is usual for such stars in the field). 
This effect decreases only weakly with distance from the BH because the increased tidal coupling 
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at lower collision velocities largely compensates for the decrease in the collision rate at lower stellar 
densities. We estimated also the stellar survival probability against head-on collisions over 10 Gyr, 
and found that it was significantly large even very close to the center. The fact that the spin-up 
is roughly constant over the volume of the inner GC implies that no large error was introduced 
by neglecting the fact that stars on non-circular orbits sample a varying stellar density over their 
lifetime. The spin-up of giant stars over the giant phase is much smaller, or order 1-2% of the 
centrifugal break-up velocity (doubling the rotation that is usual for such stars in the field). The 
effect is smaller because of their short lifetime and because of the large ratio between the collision 
velocity and the stellar escape velocity, which decreases the tidal coupling. 

These results suggest that stochastic spin-up is an important stellar effect in BH cusps. 
However, several caveats apply. This results depend on the validity of our extrapolation to the 
non- linear regime. A real star is not well represented by an ideal gas n = 1.5 polytrope, and more 
realistic and extensive hydro dynamical simulations, including the effects of radiation and nuclear 
burning, will be needed to verify our results. We note however that the n = 1.5 polytrope model is 
significantly less bound than the realistic stellar models (Table |l|). This implies that our estimates 
for the minimal peri-separation and the tidal radius may be overly conservative. In addition, 
because a smaller binding energy, larger tidal coupling coefficients (Fig. |]) and a larger moment 
of inertia go hand in hand, the spin-up (AQ ~ Ti/I) is less sensitive to the details of the stellar 
structure than any of these quantities separately. Another uncertainty in applying our results to 
any specific system, such as the GC, lies in modeling the PMF of the impactors. Generally, mass 
segregation will work towards increasing the spin-up effect by pushing the ineffective low mass 
projectiles out of the central region. Dynamical models of the evolution of the galactic nucleus 
are required to put this on a quantitative footing. We considered the survival probability of the 
stars against head-on collisions, but there are also other effects that may compete against the 
spin-up, which were not included in our estimate. Magnetic breaking, although slow, may be 
effective over 10 Gyr. Unfortunately, at present the details of this process can not be modeled 
with any certainty. The cumulative mass-loss in the course of many non-linear collisions was also 
not considered. 

Detailed predictions for the observational consequences of high rotation are outside the 
scope of this work. We limit our comments on this matter to noting that rotation lowers the 
effective temperature and luminosity of a star (Kippenhahn, Meyer- Hofmeister &; Thomas 1970| ), 



but does not significantly affect the spectral classification of a star until it is close to break-up 
(Gray |1992 ). The long term effects of rotation on stellar evolution may be more significant. We 



calculated the effective angular velocity, assuming solid body rotation. The actual distribution of 
angular momentum could become stratified over time and lead to rotational support of the core 
(VandenBerg et al. |1998[ ) or to the replenishment the hydrogen in the core by large scale deep 
mixing (Sweigart 1997] ). These effects will manifest themselves in the giant phase of the stars. 



Finally, We note that even in the GC only the giant stars can be presently observed. Rotational 
broadening in the giant spectra may be marginally detectable with high resolution spectroscopy, 
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and could bolster the case for the existence of an underlying very dense population of faint stars. 

To summarize, we have shown that MS stars in a substantial volume of the dense cusps 
around massive black holes are likely to rotate at a significant fraction of the centrifugal break-up 
velocity due to stochastic spin-up by hyperbolic tidal encounters. 



A. Linear tidal coupling coefficients for hyperbolic orbits 



Following the formalism of Press & Teukolsky (1977), the linear tidal coupling coefficients are 
expressed as 



T l (r ] )=27r 2 Y / \Qm\ 2 E \ R m 

n m=-l 



2 

m | j 



where Q n \ are overlap integrals that depend only on the stellar structure. The orbit enters in the 
term 



K, 



Wlr. 



1/2 



nlm 



TT \ 1 + fh 



:i+e) 



-1+1/2 



d(f){l + e cos i 



il-i 



cos{uj n it{4>) + m(f 



(Al) 



where iv n i is the frequency of the mode, (f) is the angular position of the impactor in a coordinate 
system centered on the target star (<j) = at periastron), and where 



'mnx arccos 1 - 



and 



1/2 



The time along the hyperbolic orbit as function of the angle <j) is given by 



(A2) 



(A3) 



-3/2 



sm< 



(1 + ef' 2 

(l + m) 1 /2 ll + ecos0e 2 -l 



e 2 -l 



-3/2 



log 



e + cos (f) + (e 2 - l) 1 / 2 sin( 
1 + e cos d> 



(A4) 



Tables |] and ^ list the mode frequencies and the overlap integrals for the solar and giant 
models that we investigate in this work. The classification of the giant modes is complicated by 
the large value of the Brunt- Vaisala frequency in the core, compared to the /-mode frequencies, 
thereby giving rise to a mix of p and g mode behavior. The corresponding values for ideal gas 
n = 1.5, 2, and 3 polytropes are given in Lee & Ostriker (1986). 



We are grateful to P. Demarque for providing us with the stellar structure model of the aUMa 
giant. 
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Type 


M 


i? 


v e 


v* 


v ohs 


K T 


I 


E b 




M Q 


R Q 


kms -1 


kms" 1 


kms -1 


mag yr 






K2(V) 


0.76 


0.75 


621 


440 


~ 2 


~ 22 > 10 10 


0.07 


1.65 


G5(III) 


2.4 


8 


340 


240 


~ 5 


~ 16 1.6 x 10 s 


0.13 


6.11 


n = 1.5 














0.21 


0.85 



Table 1: The parameters of representative MS and giant stars (Zombeck |1990 ; Lang 199 1| ). The 



mean observed equatorial circular velocity is from Gray ( |1992j ). Over the time T, the star 
is as bright or brighter than the quoted apparent K magnitude (for stars in the GC, A = 14. 5 m , 
Ak = 3.5 m ) and has a radius as large or larger than the quoted value (based on the twice solar 
metallicity stellar tracks of Schaerer et al. 1993| ). The moment of inertia /, and the binding energy 



Ej) were estimated from stellar structure models of the Sun (Christensen-Dalsgaard et al. |1996[) 
and the aUMa giant (M = 4.25 M , R = 27.4i? ) (Guenther et al. |200C| ). For comparison, 7 and 
Eb for an ideal gas n = 1.5 polytrope are also listed. 



Initial mass range 


Remnant mass 


Mass fraction 


M 


M & 




0.8-1.5 


0.6 


0.03 


1.5-2.5 


0.7 


0.08 


2.5-8 


1.1 


0.12 


8-30 


1.4 


0.03 


>30 


10 


0.01 



Table 2: Stellar remnant mass as function of initial stellar mass (Meylan & Mayor 1991 ; Timmes, 
Woosley & Weaver [1996D and its mass fraction in the continuous star forming stellar population 
model for the GC (Alexander & Sternberg 1999 ). 
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Type 


Mass range 


a 


Mass fraction 




M 






MS 


0.4-4.0 


2.35 


0.73 


WD 


0.7-1.1 


0.0 


0.23 


NS 


1.2-1.5 


2.35 


0.03 


BH 


7.0-12.0 


2.35 


0.01 



Table 3: Model for the impactor mass function in the inner GC. A df /dm am a mass distribution 
is assumed within the mass range. 



12 3 4 5 



Mode 


w 


2 
nl 


\Qnl\ 


w 


2 
nl 


\Qnl\ 


w 


2 
nl 


\Qnl\ 


w 


2 
nl 


\Qnl\ 


P6 


1.28 


:+2) 


2.38( 


-2) 


1.39 


:+2) 


1.12( 


-2) 


1.49 


:+2) 


5.39( 


-3) 


1.59 


:+2) 


2.74( 


-3) 


P5 


9.70 




3.26( 


-2) 


1.06 


:+2) 


1.58( 


-2) 


1.14 


:+2) 


8.00( 


-3) 


1.23 


:+2) 


4.28( 


-3) 


Pi 


7.02 




4.67( 


-2) 


7.76 




2.34( 


-2) 


8.48 


:+i) 


1.22( 


-2) 


9.18 




6.82( 


-3) 


P3 


4.79 




6.85( 


-2) 


5.36 




3.48( 


-2) 


5.90 




1.91( 


-2) 


6.45 




1.14( 


-2) 


P'2 


2.95 


:+i) 


1.06( 


-1) 


3.33 




5.61( 


-2) 


3.72 




3.32( 


-2) 


4.12 




2.16( 


-2) 


Pi 


1.70 




1.471 


-1) 


1.82 




9.18( 


-2) 


2.01 




6.36( 


-2) 


2.23 


:+i) 


4.66( 


-2) 


f 


1.39 




1.39( 


-1) 


1.65 




5.54( 


-2) 


1.77 


:+i) 


8.62( 


-3) 


1.84 


:+i) 


7.56( 


-4) 


91 


9.12 


;+o) 


4.58( 


-2) 


1.17 




8.47( 


-3) 


1.35 




4.80( 


-3) 


1.47 




2.40( 


-3) 


92 


6.92 


:+o) 


1.13( 


-1) 


8.94 


:+o) 


4.46( 


-2) 


1.08 




1.64( 


-2) 


1.23 




7.29( 


-3) 


93 


5.11 


:+o) 


1.07( 


-1) 


7.02 


:+o) 


8.99( 


-2) 


8.58 


:+o) 


3.81( 


-2) 


1.00 




1.56( 


-2) 


94 


3.84 


;+o) 


8.31( 


-2) 


5.81 


;+o) 


1.10( 


-1) 


7.13 


;+o) 


1.00( 


-1) 


8.35 


;+o) 


4.85( 


-2) 


95 


2.92 


;+o) 


5.66( 


-2) 


4.75 


;+o) 


7.20( 


-2) 


6.35 


:+o) 


1.04( 


-1) 


7.51 


;+o) 


1.29( 


-1) 


96 


2.31 


;+o) 


3.91( 


-2) 


3.84 


;+o) 


4.36( 


-2) 


5.38 


;+o) 


4.76( 


-2) 


6.76 


;+o) 


5.20( 


-2) 


97 


1.82 


:+o) 


2.771 


-2) 


3.13 


:+o) 


2.90( 


-2) 


4.45 


:+o) 


2.73( 


-2) 


5.76 


:+o) 


2.31( 


-2) 


9s 


1.49 


;+o) 


1.99( 


-2) 


2.59 


;+o) 


2.04( 


-2) 


3.76 


;+o) 


1.84( 


-2) 


4.93 


;+o) 


1.48( 


-2) 


99 


1.23 


;+o) 


1.44( 


-2) 


2.19 


;+o) 


1.47( 


-2) 


3.20 


;+o) 


1.30( 


-2) 


4.24 


:+o) 


1.03( 


-2) 


9w 


1.04 


;+o) 


1.07( 


-2) 


1.85 


;+o) 


1.12( 


-2) 


2.76 


;+o) 


1.02( 


-2) 


3.69 


;+o) 


8.20( 


-3) 



Table 4: Mode frequencies and overlap integrals for the solar model of Christensen-Dalsgaard et 
al. ( gggg ). 
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e = 1 





Fig. 1. — The first few orders of 7} for the Solar stellar structure model, the giant model and the 
ideal gas n = 1.5 polytrope. Top: parabolic (e = 1) encounters. Bottom: hyperbolic (e = 10) 
encounters. 
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Fig. 2. — The energy transferee! from the orbit to the star in a parabolic collision between a point 
mass and an ideal gas n = 1.5 polytrope of equal mass (m = 1, E Q = 0) as function of the peri- 
separation r p . The SPH results from this work are compared with the higher resolution SPH results 
of Rasio & Shapiro ( |1991| ) and with the predictions of the linear theory (Eq. |12| with Cnl = !)■ 



I 




2 






3 






4 




Mode 


W 


2 
nl 


\Qnl\ 


w 


2 
nl 


\Qnl\ 


w 


2 
nl 


\Qnl\ 


P8 


4.28 


:+i) 


1.25(-2) 


4.53 




1.62(-2) 


4.83 




1.86(-2) 


P7 


3.66 


:+i) 


1.66(-2) 


3.40 


:+i) 


2.82(-2) 


3.67 


:+i) 


3.18(-2) 


P6 


3.08 




2.29(-2) 


2.21 




6.09(-2) 


3.15 


:+i) 


4.12(-2) 


P5 


2.19 




4.35(-2) 


1.60 




1.02(-1) 


2.62 




5.86(-2) 


Pi 


1.45 




8.40(-2) 


1.11 




1.66(-1) 


1.68 




1.13(-1) 


P3 


1.27 


:+i) 


1.05(-1) 


9.24 


:+o) 


2.47(-l) 


1.01 


:+i) 


2.27(-l) 


P2 


8.12 


:+o) 


2.36(-l) 


6.55 


;+o) 


1.26(-1) 


6.81 


;+o) 


3.59(-2) 


Pi 


6.10 


:+o) 


2.68(-l) 


4.45 


:+o) 


9.06(-2) 


5.02 


:+o) 


3.42(-2) 


f 


3.57 


:+o) 


1.84(-1) 


3.39 


;+o) 


8.86(-2) 


2.19 


;+o) 


3.69(-2) 



Table 5: Mode frequencies and overlap integrals for the giant model of Guenther et al. (|2000 ). 
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X— Y projection 



— | t 2.84E+00 

- step 450 

- M*1 1.000 

- M*2 1.000 

- Mub 0.000 



t 3.75E+00 
step 650 

M»1 0.925 
M*2 0.952 
Mub 0.146 



-2 2 

X— Y projection 



— | t 5.28E+00 

- step 550 

- M*1 0.875 

- M*2 0.B77 

- Mub 0.250 



t 5.68E+00 
step 900 

M*1 0.944 
M*2 0.945 
Mub 0.112 



Fig. 3. — A sequence of snapshots from an SPH simulation of a deep collision between two equal 
mass stars (modeled as ideal gas n = 1.5 polytropes with N = 1024 particles each). The orbital 
parameters of the encounter are r p = 1.0 and E Q = 2. At t = 20 after the periastron passage (not 
shown here), the stars have lost 5% of their mass each, were spun-up by = 0.025 and acquired 
a moment of inertia / = 0.47, more than twice the initial value. 
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Fig. 4. — The response of a star [fa = 1, modeled as an ideal gas n = 1.5 polytrope) to deep 
non-linear hyperbolic collisions (E Q = 2) with another star, for both the case where the impactor 
(modeled as a point mass) is compact relative to the target star (circles) and the case where 
the impactor (modeled as an ideal gas n = 1.5 polytrope) is of the same size as the target star 
(triangles). The SPH results (thin lines with symbols) are compared with the analytic calculations 
without the non-linear correction (thick lines) for AL (top panel) and AO, (middle panel). Also 
shown in the bottom panel are 1 — AM*, the mass that remains bound to the star after the collision 
(thin lines with filled symbols) and /, the star's moment of inertia after the collision (thin lines 
with open symbols). 
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X-Y projection 
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X-Y projection 



t 6.0SE+00 
step 3200 
d 10.339 
M* 0.811 
Mub 0.033 
Mpm 10.156 




Fig. 5. — A sequence of snapshots from an SPH simulation of a collision at the tidal limit (Eq. |l6|) 
between a in = 10 BH (modeled as a point mass) and a MS star (modeled as an ideal gas n = 1.5 
polytrope). The dots represent SPH gas particles that remain bound to the star, and the circles 
represent those that are lost. The line traces the trajectory of the BH in a frame where the target 
was initially at rest. The orbital parameters of the encounter are f p = ft = 2.0 and E Q = 1.0. In 
spite of the significant tidal stretching of the star (bottom left panel), the star ultimately relaxes 
after losing 25% of its mass. At t = 124 (bottom right panel) the stellar mass that has settled back 
to within f = 2 of the stellar core rotates at 20% of the centrifugal break-up velocity. The moment 
of inertia of the mass within f = 2 is I = 0.45, more than twice its initial value. 
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Fig. 6. — A comparison of the DFs of the relative stellar velocity (w = V2/^ 1 ^ 2 ) in a Maxwellian 
DF and of a e p DF of relaxed low- mass stars very near a black hole (Bahcall & Wolf 1977). Both 
DFs have the same ID velocity dispersion and the same mean orbital energy. 
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Fig. 7. — The mean stellar spin Af2 of a MS star as function of distance from the BH in the GC 
assuming a Salpeter PMF and the remnant mass fractions listed in table 0. The survival probability 
against a destructive collision with f p < 0.5 over 10 Gyr (Eq. |i~9| ) is also plotted for a mean impactor 
mass of (m) = 0.9 Mq. 



